Parametric excitation of a Bose-Einstein condensate in a ID optical lattice 
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We study the response of a Bose-Einstein condensate to a periodic modulation of the depth of an 
optical lattice. Using Gross-Pitaevskii theory, we show that a modulation at frequency fi drives the 
parametric excitation of Bogoliubov modes with frequency fi/2. The ensuing nonlinear dynamics 
leads to a rapid broadening of the momentum distribution and a consequent large increase of the 
condensate size after free expansion. We show that this process does not require the presence of a 
large condensate depletion. Our results reproduce the main features of the spectrum measured in 
the superffuid phase by Stoferle et at, Phys. Rev. Lett. 92, 130403 (2004). 
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Exploring the effects of interaction and low dimension- 
ality is one of the central themes in current research on 
ultracold Bose gases. Optical lattices are ideal tools in 
this sense. A major achievement was the demonstra- 
tion that, by increasing the depth of an optical lattice, a 
quantum phase transition occurs from the superfluid to 
a Mott insulating phase 0. More recently, Stoferle et 
al. 0, Q made an effort to characterize this phenomenon 
by measuring the excitation spectrum. The atoms are 
loaded into the ground state of an harmonic trap plus 
an optical lattice, and they are subsequently excited by 
modulating the lattice depth. Then the gas is released 
from the trap and the width of the expanding cloud is 
taken as a measure of the excitation energy. The quan- 
tum phase transition shows up as the crossover from a 
broad continuum of excitations in the superfluid phase 
to a discrete spectrum in the Mott insulating phase. 

The broad resonance in the superfluid regime was un- 
expected. In fact, the Gross-Pitaevskii (GP) theory pre- 
dicts that, in the linear response regime, a modulation of 
the lattice depth cannot excite the gas 0. For this rea- 
son, it was conjectured that the measured spectra might 
be the effect of strong interactions, not included in the 
GP theory, which yield a significant quantum depletion of 
the condensate Here we show that the observed res- 

onance is caused by the parametric amplification of Bo- 
goliubov states and the subsequent nonlinear dynamics, 
which leads to the broadening of the momentum distribu- 
tion. We describe both processes by using the GP theory 
and assuming the initial condensate to be out of equilib- 
rium. A tiny initial occupation of Bogoliubov modes is 
found to be enough to produce a strong response at high 
frequency, while larger fluctuations are needed to explain 
the observed spectrum at low frequency. 

We simulate the dynamics of the condensate by solving 
the GP equation 6] for an axially symmetric condensate 
of 87 Rb, similar to typical "tubes" in the experiments 
of Ref . 2] . Since the transverse and axial trapping fre- 
quencies are such that uj±/lu z 3> I, one can factorize the 
order parameter in the product of a Gaussian radial com- 
ponent of z-and i-dependent width, a{z, t), and an axial 



wave function ^(z,t). The GP equation yields 
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and a = a ± (l + 2a7V|*| 2 ) 1 / 4 . Here a x = [h/ '(rawi)] 1 ' 2 , 
N is the number of particles, g = Anh a/m, with a the 
s-wave scattering length and m the atomic mass. The 
axial potential V is the sum of the harmonic trap and 
the optical lattice: 

V{z,t) = { Jn /2)uj 2 z z 2 +sE R [l+Aam(nt)}sm 2 {qBz) , (2) 

where q R = n/d is the Bragg wave vector, d is the lattice 
spacing and s is the lattice depth in units of the recoil 
energy E R = Ti 2 q 2 B /2m. To simulate a single tube of 
Ref. [3| , we use d = 413 nm, uj z = 2ir x 84.6 Hz and 
lu ± = 2n x 36.5 KHz, yielding E R /h = 2tt x 3.34 KHz. 
We first calculate the ground state for given N and s, by 
solving the stationary GP equation with A = 0. Then we 
solve the time-dependent equation with A ^ for a 
modulation time i m . Finally, we switch off the external 
potential and we let the condensate expand for 25 ms || . 

In Figs. 1 and 2 we report the results for N = 100, 
s = 4, A = 0.2 and Q = 1.33E R /h. The solid line in 
the upper panel of Fig. shows the width of the density 
distribution at the end of the expansion as a function of 
the in-trap modulation time t m 9J. Before about 20 ms 
nothing seems to occur. Then the width starts increas- 
ing and eventually reaches values of about 10 times the 
size of the expanding ground state. The density and mo- 
mentum distributions of the in-trap condensate at the 
instants A, B, C, and D are plotted in Fig. 2. At the 
early time A, both are still indistinguishable from those 
of the ground state. In contrast, the momentum distri- 
bution in B clearly features new sharp peaks at ±g, with 
q = 0.56gs and the density is modulated on a lengthscale 
of about 4 sites, accordingly. By the time C, the popula- 
tions at ±q have grown large. Besides, the lateral peaks 
at ±2qs, associated with the density modulation of the 
ground state in the lattice, are smaller and contributions 
from momenta throughout the first and second Brillouin 
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zones start showing up. At D, the system hardly bears 
any resemblance with the ground state and its broad mo- 
mentum distribution explains the increase of the width 
of the expanding cloud shown in Fig. ^ 
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of parametric instabilities (see e.g., [llj and references 
therein) . In our case, the modulation of the lattice depth 
at frequency fl causes the instability of modes with fre- 
quency w(q) = Q/2. To demonstrate this, we perform 
simulations at various f2 and, from each one, we extract 
the positions ±q of the first extra peaks which grow in the 
momentum distribution. Given these values of q, in Fig. 3 
we plot the frequencies fl/2 as a function of q (points), 
together with the dispersion law of the lowest Bogoliubov 
band, ui(q) (solid line). The latter is calculated by keep- 
ing the same optical lattice and transverse trap, but ne- 
glecting the weak axial confinement; the number of atoms 
per lattice site is taken to be equal to the average num- 
ber of atoms per site in the actual trapped condensate 
[lil Haj ] . The agreement between points and solid line 
clearly indicates that the dynamics of the condensate is 
directly related to the spectrum of Bogoliubov modes. 



FIG. 1: Upper panel: the width of the density distribution 
<2 2 ) 1/2 after 25 ms of free expansion is plotted as a function of 
the modulation time t m , as obtained with £1 — 1.33Er/?i and 
A = 0.2; Lower panel: population N q of the parametrically 
excited mode with q = 0.56gs and uj = fi/2. The different 
lines refer to different initial states. Solid line: ground state; 
dashed: ground state plus a weak Bragg pulse; dotted: ground 
state plus a small white noise. 
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FIG. 2: Density (left) and momentum (right) distributions at 
different times A, B, C and D given in Fig. 

This behavior is a consequence of the parametric in- 
stability of Bogoliubov modes. A parametric instability 
corresponds to the exponential growth of certain modes 
of a syst em induced by the periodic variation of a param- 
eter [lOj. Systems governed by a nonlinear Schrodinger 
equation, as the GP equation, can exhibit a wide variety 
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FIG. 3: For each modulation of given frequency fi, we draw 
a point at (q,co), where q is the wavevector of the parametri- 
cally excited mode and ui = Q/2. The solid line is the lowest 
Bogoliubov band. 

In order to be parametrically amplified, the modes at 
±g must preexist at t = 0. This means that our initial 
state is not a pure ground state. Indeed some numerical 
noise is always present and yields a very small occupa- 
tion of Bogoliubov modes. To understand this point, we 
define the fraction of atoms, N q (t), which contributes to 
the peaks at ±q Q and plot it in Fig. 1. The solid 
line starts at about 5 x 10~ 6 which, for the parameters 
of this simulation, coincides with the iq components of 
the exact ground state. These nonvanishing components 
are due to the finite axial size of the condensate; they 
do not grow in time and must not be confused with the 
population of the Bogoliubov modes. The seed popula- 
tion is instead provided by the numerical noise and is 
much less than 10~ 6 at t = 0. Then it grows exponen- 
tially and it becomes visible in Fig. 1 after about 15 ms, 
when it exceeds the ground state level. It eventually sat- 
urates at a value of the order of 10 _1 . The time scale 
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of this growth is much longer than the lattice modula- 
tion period. The growth rate that we extract from our 
simulations is linear in il and A within our accuracy. 
This linear dependence agrees with the results of a semi- 
analytic model which describes the mean-field coupling 
between the ground state and ±o excitations within a 
quasiparticle projection method [15|. 

The time needed to saturate depends on the amount 
of seed. To explore this dependence, we mimic an initial 
depletion by introducing some extra "noise" in a control- 
lable way. For instance, just before t — we excite the 
condensate with a short and weak Bragg pulse in reso- 
nance with the Bogoliubov mode at q. Alternatively, we 
simply add random noise with uniform (white noise) or 
thermal-like distributions. The effects of the Bragg pulse 
and white noise are shown in Fig.^ where the seed is cho- 
sen to be ~ 10~ 5 in both cases. The main effect of this 
extra seed is to anticipate the instant at which the para- 
metric excitation becomes noticeable. Yet, the growth 
rate and the level at which (z 2 ) 1 / 2 and N q saturate re- 
main the same as for the solid line, which corresponds 
to less than 10~ 15 seed excitations. This means that, for 
this f2, the amount of seed which is enough to trigger the 
parametric instability is indeed very small. In the ex- 
periments the system is certainly not in its groundstate 
at t = due to thermal and quantum depletion, and to 
excitations produced when loading the condensate into 
the lattice 16]. Our calculations show that even a very 
tiny deviation from the groundstate can produce a strong 
response on the timescale of the experiment. 

The evolution leading up to the results in Figs. 1 and 
2 not only involves the parametric excitation of the ±q 
modes, but also the subsequent broadening of the mo- 
mentum distribution, as in C and D. This is due to the 
nonlinear mean-field dynamics of the condensate gov- 
erned by the GP equation. This nonlinear dynamics is 
triggered by the growth of the parametrically unstable 
modes, but its outcome is little affected by changing 
or A. Even more, it is hardly influenced by the external 
perturbation; in fact, switching off the lattice modulation 
at time C has very little effect on the outcome at time D 
or later. This implies that the response of the condensate 
as a function of £1 mainly reflects the il-dependence of the 
growth rates of the parametrically excited modes. More- 
over, the time scale for the nonlinear broadening of the 
momentum distribution is also relatively long, so that the 
increase of the width (z 2 ) 1 / 2 occurs on a scale of several 
tens of ms, as in Fig. 1. This agrees with the behavior 
observed in the experiments (see Fig. 3 of Ref. [r^h 

As seen in Fig. 3, the main effect of changing the lattice 
modulation frequency is that of picking out a different 
Bogoliubov mode for parametric excitation. In the limit 
£1 — » the growth rate vanishes and the width (z 2 ) 1 / 2 at 
a fixed t m becomes equal to the width of a condensate 
which expands without being modulated. The same lim- 
iting value is obtained on the opposite side, namely at £1 



close to twice the Bogoliubov bandwidth, 2w(<7b), where 
parametric excitations are not possible due to the pres- 
ence of the gap in the Bogoliubov spectrum at the zone 
boundary. Therefore the response of the system takes 
the form of a broad peak in the range < £1 < 2cj(<j , b). 
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FIG. 4: Width <z 2 ) 1/2 of the expanded array of condensates 
as function of fi. Top panel: experimental data of Ref. 0, 
for s — 4; Solid line in bottom panel: our simulations starting 
from a small white noise. Empty circles and solid squares: 
same but with a much larger white noise and thermal-like 
noise, respectively (see text for explanation). 

In order to compare our results with the experimental 
observations of Ref. 0, a simulation with a single tube 
is not enough. The actual condensate is made of an ar- 
ray of tubes having different population, N, and hence 
different Bogoliubov bandwidths. Given the total num- 
ber of atoms in the array and the trapping frequencies, 
one can estimate the population of each tube by calcu- 
lating the "average" density distribution within a local 
density approximation as in Refs. 003. Then we per- 
form several simulations for different values of N and we 
sum their contributions to (z 2 ) 1 / 2 using the "average" 
density distribution as a weighting function. The central 
condensate is taken to have N = 90 which is expected 
to be reasonably close to reality. The results are shown 
in the lower panel of Fig. 4. The solid line is the result 
of GP simulations with a white noise as a seed excita- 
tion, corresponding to N q ~ 10~ 3 at t = 0. This small 
seed is enough to produce a broad response which quali- 
tatively agrees with the experimental one (upper panel). 
The height of the peak is of the same order and its slope 
at high frequency is very similar. Within our approach, 
this slope originates from the fact that the response van- 
ishes at twice the Bogoliubov bandwidth, 2a;(gs), but 
this bandwidth is different for tubes with different N in 
the array and this provides a smoothing of the high fre- 
quency tail. By repeating the calculation for a deeper 
lattice (s = 6, instead of 4), we find that this tail shifts 
to lower frequencies by about 10-15%. This agrees with 
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the experimental observations and is consistent with the 
calculated decrease of the bandwidth. 

On the other side of the peak, at low frequency, the re- 
sponse turns out to be more sensitive to the seed. In or- 
der to reproduce the observed response around 2000 Hz, 
we have to start with a condensate in which a significant 
number of atoms populate the lowest Bogoliubov modes. 
This also yields an increase of the axial width of the ex- 
panded condensate without modulation. The effect of a 
larger seed is represented in Fig. 4 by empty circles and 
solid squares, which are the results of calculations with 
a white noise and a thermal-like seed, respectively. The 
thermal seed is simulated by random fluctuations, where 
the average occupation number of Bogoliubov excitations 
follows a Bose distribution at a given temperature. The 
dispersion u>(q) of a uniform condensate, calculated as 
for Fig. 3, and the corresponding quasiparticle ampli- 
tudes are used to relate the occupation of the quasipar- 
ticle states to the momentum distribution of the atoms. 
Both the temperature in the thermal-like seed and the 
level of the white noise are chosen in such a way that the 
width without modulation is close to the experimental 
value, (z 2 ) 1 / 2 ~ 25/im. This implies a total depletion 
of the order of 20% in the case of white noise and the 
temperature fc^T ~2x 10^ 3 huj± for thermal noise. The 
latter value is compatible with the experimental estimate 
ksT < 6 x 10~ 3 ?kj^ given in Ref. [19(. As one can see, 
in both cases the overall response is higher and broader, 
and the qualitative agreement with the experiments is 
improved. This result is also compatible with the ob- 
servation made in |3j that the excitation efficiency de- 
creases when the transverse confinement is reduced and, 
hence, the initial depletion of the condensate is reduced 
as well. A quantitative comparison is however difficult 
at this level. In the experiment a thermalization process 
is performed after the modulation and before the expan- 
sion, by lowering the depth of the transverse optical lat- 
tice and letting the tubes interact for a while. This is 
not included in our simulations, since it would require a 
full 3D calculation. Also the procedure of averaging over 
many tubes, some of them with only a few atoms, limits 
the accuracy of the comparison. Nevertheless, our anal- 
ysis suggests that a characterization of classical and/or 
quantum fluctuations can indeed be possible by appro- 
priately choosing the trapping geometry and the exper- 
imental procedure and using the parametric resonances 
for a selective amplification of initial fluctuations [lfij . 

In conclusion we have shown that the broad resonance 
observed in condensates subject to a periodic modula- 
tion of the optical lattice can be interpreted as the effect 
of a parametric instability of Bogoliubov modes. The 
occurrence of this type of instability in Bose-Einstein 
condensates is an interesting property of these systems, 
which originates from their superfluid nature and the ex- 
istence of undamped collective excitations. The paramet- 
ric instability and the subsequent nonlinear dynamics are 



nicely accounted for by the GP equation. Since our cal- 
culations are purely mean-field in nature, they are not 
conclusive in excluding possible many-body processes in 
the explanation of the experimental observations. Ex- 
ploiting all predictions of GP theory is however essential 
in view of investigations on beyond mean-field effects due 
to strong interactions and low dimensionality. 
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